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ABSTRACT 



Super-hot looptop sources, detected in some large solar flares, are compact 
ms> \ sources of HXR emission with spectra matching thermal electron populations ex- 

C/^ ! ceeding 30 megakelvins. High observed emission measure, as well as inference of 

^ ! electron thermalization within the small source region, both provide evidence of 

I ! high densities at the looptop; typically more than an order of magnitude above 

^ ! ambient. Where some investigators have suggested such density enhancement 

c/5 , results from a rapid enhancement in the magnetic field strength, we propose an 

alternative model, based on Petschek reconnection, whereby looptop plasma is 
heated and compressed by slow magnetosonic shocks generated self-consistently 
through flux retraction following reconnection. Under steady conditions such 
Tj- \ shocks can enhance density by no more than a factor of four. These steady shock 

relations (Rankine-Hugoniot relations) turn out to be inapplicable to Petschek's 
model owing to transient effects of thermal conduction. The actual density en- 
O \ hancement can in fact exceed a factor of ten over the entire reconnection outflow. 

An ensemble of flux tubes retracting following reconnection at an ensemble of dis- 
tinct sites will have a collective emission measure proportional to the rate of flux 
'k/. \ tube production. This rate, distinct from the local reconnection rate within a 

H \ single tube, can be measured separately through flare ribbon motion. Typical 

flux transfer rates and loop parameters yield emission measures comparable to 
those observed in super-hot sources. 

Subject headings: MHD — Sun: flares — Sun: magnetic fields 



1. Introduction 

Solar flares are best known for releasing large amounts of magnetic energy. They are 
equally remarkable, however, for large increases in plasma density accompanying the energy 
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release. During a large flare the coronal density can increase by two or three orders of 
magnitude, from an ambient value like Ue ~ lO^cm"^, to rig ~ 10^^ or 10^^ cm~^ (see eg. 
Moore et al.lll980l ) . The temperature of this dense plasma is typically one order of magnitude 
larger than ambient conditions (Tg ~ 20 MK) so the magnetic field must confine plasma 
pressure three to four orders of magnitude above normal. The coronal magnetic field is strong 
enough to confine this high-pressure plasma in the two directions perpendicular to field lines, 
but the plasma is unconfined along the field lines. It is commonly believed that transport of 
mass and energy along the field lines drives a process known as chro mospheric evaporation 
to su pply the extra mass responsible for the density enhancement (jAntiochos fc Sturrock 
19781 ). Energy is added to the chromosphere causing it to expand into the corona, thereby 
diluting the coronal energy over more mass, lowering its temperature but raising its density. 

In addition to hot loops, filled by chromospheric evaporation, some flares include an 
over-dense "knot " of plasma at the loop apex known as looptop source (LTS, see review 



by iFletcheii Il999l ) . It is possible that many m ore flares have LTSs t hat cannot be imaged 
by current methods limited in dynamic range (jPetrosian et al.ll2002l ). An LTS is generally 
assigned a high plasma density due to its radiative output in hard X-rays (HXRs) by which 
the sources are most commonly imaged. Its locati on atop the IootJJ s uggests the LTS is 
a direct manifestation of magnetic energy release (JMasuda et al.lll994l ). If they are thus 
assumed to arise early in the energy release sequence, it would seem problematic to invoke 
evaporation as the source of excess density in LTSs. M oreover, the source appears physically 
separated from the chromosphere (I Jiang et al.ll2006l ). suggesting its mass is not supplied 
from there. 

A subset of LTSs appear to be dense enough that their HXR emission arises from an elec- 
tron population which is convincingly thermal (i.e. Maxwellian). The tempe ratures of these 
super-hot (SH) sources typically exceed 30 M K and sometimes reach 50 MK (JLin et al.lll981 



19851 : iKosugi et al.l 1 1994 : iNitta fc Yajil 119971 ) . One indication of their density comes from 



the emission measure (EM) of the thermal brehmsstrahlung, typically exceeding 10 cm . 
Attributing a large filling f actor (often unity) t o the apparent source provides density y al- 



ues as high as rig ~ 10 fIVeronig et al.l l2006l : iLongcope et al.l l2010l : ICaspi fc Lid 12010 



smaller filling factor would demand an even higher density). Independent evidence for these 
high densities comes from the inference of a Maxw ellian electron distribution, supported 



at high confidence by high resolution HXR spectra (JLin et al.l Il98ll : iLongcope et al.l 12010 



Caspi fc Linll2010l ). Electrons must collide with one another numerous times to fully ther- 



^Some investigators prefer the term "above-the- looptop" in recognition of its separation from loops visible 
in soft X-ray or EUV images. In all likelihood the LTS is at the top o/some magnetic loop, even if that loop 
is not (yet) visible at longer wavelengths. 
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malize. Accommodating this within a source extending only 10 Mm requires electrons with 
mean free paths ^e ^ 1 Mm along the field line. Electrons at a typical SH temperature of 
Te = 35 MK have a collisional mean free path, 



100 Mm 



Tie, 



109cm-3 



-1 



along the magnetic field line, so thermalization within the LTS demands ng > 10 cm . 

Generally, LTSs show a non-thermal electron p opulation (e.g. a power-law ), similar to 
footpoint sources, during a flare's impulsive phase ([Alexander &: Metcalil 119971 ) . Superhot 
thermal sources typically appear after the ir npulsive phase and persist ~ 10 minutes or more, 
far longer than the conductive cooling time (j Jiang et al.ll2006l ). Some HXR spectral analyses 
show the SH source co-existing with a secon d thermal population at Tp ^ 20 M K and much 



higher emission measure (5-20 times higher. iLongcope et al.ll2010l : ICaspi fc Linll2010[ ). This 
lower-temperature component closely matches the soft X-ray light-curve, from GOES for 
example, commonly attributed to chromospheric evaporation. This re-affirms a role for the 
SH component as a fiare stage prior to evaporation, thus demanding a different mechanism 
for enhancing its density. 

To date, little effort has been devoted theoretically explaining the high densities in LTSs. 
Non-thermal sources are sometimes a ssumed to be filled by chr omospheric evaporation driven 
by impact of precipitating particles (jVeronig fc Brownll2004l ). Doing so tacitly assumes an 
energy release process able to persist long enough on a single field line to interact with 
its own effects: chromospheric evaporation. This seems at odds with the favored energy 
release mechanism, magnetic reconnection, whereby a single field line will reconnect once in 
a single instant. Resolution of this puzzle seems to require a self-consistently coupled model 
of particle acceleration and magnetic reconnection; such a model is still still being pursued. 

Since SH looptop sources are thermal, and show little evidence for non-thermal par- 
ticles, they might be understood without an understanding of p article ac c elerat ion. Most 
widely accepted models of magnetic reconnection, such as that of lPetschekl (119641 ). use fiuid 
equations and are therefore directly applicable to SH sources. Several LTS models proposed 
so far have invoked an increa se in magnetic fields strength, a "collapsing trap " , to enhance 
the d ensity and temperature (ISomov fc Kosugil 119971 : iKarlicky fc Kosugil 12004 : ICaspi fc Lin 
20101). To match observed densities these often require field strengths to increase by factors 
of 30-100. These models were not based on reconnection solutions, in fact, assuming a field 
strength increase would seem to contradict spontaneous reconnection scenarios where energy 
is released by reducing magnetic energy. The notable exception is the fast magnetosonic ter- 
mination shock which sometimes forms where the reconnection outfiow jet encounters an 
"obstacle" , such as the arcade of previously reconnected loops. 
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Fast magnetoso nic termination shock s have been a common feature in reconnection- 



driven flare models (IForbes fc PriestI Il983l : iForbed Il986l : iTsuneta et al.l 119971 : lAurass et al. 
2OO2I : IVrsnak fc Skenderl |2005| ) . The magnetic field strengt h increases across a sh ock of 



this type, making it suitable for Fermi particle acceleration (ITsuneta &: Naitdll998). The 
densi t y also increases across th e shock, but typically by no more than a factor of two (JForbes 



1986 



Vrsnak fc Skenderl l2005l ). This modest density enhancement (far short of the factor 



of 100 often seen in loop-top sources) coupled with the low emission measure expected in a 
small structure, makes the fast magnetosonic termination shock an unlikely candidate for a 
collapsing trap or a LTS. 

The principal element in Petschek's fast reconnection model is a slow magnetosonic 
shock, or slow shock, which heats and accelerates the plasmao Since field strength decreases 
across shocks of this type they are unsuitable for Fermi acceleration. They do, however, com- 
press the plasma considerably more than the fast magnetosoic termination shock can. Re- 
connection between anti-parallel fields generates sl ow shocks at the s witch-off limit, wherein 
the maximum possible density enhancement is 2.5 (JForbes et al.lll989l ). Reconnecting skewed 
fields (i.e. any angle other than anti-parallel) can, however, result in density enhancements 
as great as four, according to standard shock models. This compressed and heated plasma 
forms a long jet extending away from the reconnection site which may have considerably 
greater emission measure than the small fast magnetosoic termination shock. For these rea- 
sons it seems plausible for LTSs to be manifestations of the plasma heated and compressed 
by slow shocks according to Petschek's reconnection model. 



Recently iLongcope et al.l (I2OIOI ) modeled a SH-LTS observed by RHESSI using slow 
shocks from a 3-dimensional time- dependent reconnection theory. The post-shock tempera- 
ture follows directly from the angle between re connecting fields l ines, which follows in turn 
from a model of the pre-fiare magnetic field. iLongcope et al.l (120101 ) found the observed 
SH temperature T ~ 40 MK to be consistent with the angle from their pre-fiare magnetic 
model: A6 = 70° (anti-parallel fields have A6 = 180°). The time-dependent model assumed 
reconnection occurred sporadically in numerous small patches within the current sheet sepa- 
rating the skewed fields. Each patch produced a single flux tube whose subsequent retraction 



^Petschek's 1964 paper contains several pioneering elements which were extensively studied in subse- 
quent investigations under the generic term "Petschek reconnection" . It derived a steady external solu- 
tion matched to a resistive internal solution to obtain its well-kno wn reconne ction rate. This solution was 

Vasvliunaj 



improved and genera lized over following decades (jSonnerup 



1970; 



1975; 



Soward fc Priest 



1982t 



Priest fc Forbeslll986r). The solution is, ho wever, modified when the electric field is not resistive or is unsteady 
(jHevn fc Semenovlll996l : iNitta et al.ll2001 ). Petschek's work was also the first to recognize t he generation of 



slow shocks by reconnect ion, which turn out to be inevitable in most fast reconnection modes (|Semenov et al. 



19831 : lErkaev et al.ll2000l ). It is to this aspect of Petschek's work which we exclusively refer hereafter 
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through the sheet generated a plug of SH plasma between traveling slow shocks. Each plug 
emitted for ~ 8 sec before confining inflows ceased and its own pressure "disassembled" it. 
The single observed LTS was actually composed of ~ 30 separate, unresolved plugs in flux 
tubes created by separate patches across the sheet. The plug collection had a total emission 
measure proportional to the rate of patch-production which was in turn proportional to the 
mean reconnection rate. The latter was measured using the motion of chromospheric flare 
ribbons and found to be consistent with the emission measure of the observed LTS. Each 
reconnected flux tube was later observed at lower temperature as a distinct post-fiare loop 
in TRACE 171A images. The rate of loop appearances was also deemed consistent with the 
the inferred patch-production rate. 

While the compression ratio is greater for a slow shock than for the fast magnetosoic 
termination shock (4 vs. 2), it is still signific antly below the l evel i nferred for observations 
of looptop sources: 10-100. This fact forced iLongcope et al.l ( l2010l ) to assume a very high 
ambient density {rieo — 8 x 10^° cm~^) in order to match the observed emission measure. This 
same problem faces all models invoking shocks to compress plasma. The aforementioned 
compression ratios, used by all analytical reconnection models, come from standard MHD 
shock models, based on conservation laws across a steady-state jump (see iPriest fc Forbes 
2000l . for example). One possible exception occurs when radiation can cool t he plasma faster 
than compression heats it — so called radiative shocks (IXu fc Forbes! Il992l ). Unfortunately, 
the radiative loss function is so low at SH temperatures that even a density of rie = 10^^ cm~^ 
cannot produce effective cooling (xrad — 8 minutes!). Thus it would seem impossible to invoke 
a radiative shock in a SH looptop source. 

Shocks generated by magnetic reconnection are, however, so far from steady state 
conditions that standard compression ratios, so often used, are not actually apphcable. 



Therm al conduction a 
shock (Thomas! 1944 : 



ong the m a gnetic fie l d line will generate a heat fron t ahead oi the 



Grad 1951; Germain 1960; Forbes fc Malherbe 1986; Kennel 1988 



Yokoyama &: Shibata!!l997! ). Since thermal conduction does not affect net energy conserva- 
tion, the heat front does not change the steady-state density ratio. Under coronal conditions 
thermal conductivity is extremely large and the heat front would be enormous by the time 
it ach ieved steady state - far larger than any conceivable flare loop (JGuidoni fc Longcope 
2010l ). This forces the conclusion that solar flare shocks ne v er ach ieve their steady state. 
Recent simulations of these shocks by !Guidoni fc Longcope! (|2010! ) show density enhance- 
ments by factors 6-8 during the transient phase following reconnection. More careful study, 
herein presented, reveals there is no upper bound to the transient density enhancements and 
that factors of 10-100 are reasonable for a large solar flare. We therefore propose that the 
large plasma densities in LTSs are the direct result of compression in slow shocks of the fast 
reconnection liberating the energy. 
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One advantage of the model here proposed is that it directly couples the density en- 
hancement to the energy release. In models of fast magnetic reconnection, magnetic energy is 
liberated primarily by shortening field lines following topological change within a very small 
diffusion region. Except in cases of purely anti-parallel reconnection, the field strength de- 
creases only slightly and thus the release cannot be interpreted as field annihilation. The field 
line shortens at the local Alfven speed which is far larger than the hydrodynamic sound speed 
(/3 ^ 1). The plasma within the shortening flux tube is therefore compressed at super-sonic 
speeds, generating the slow shocks which heat the plasma. This is a time-dependent process 
on any given field line, so the shock behavior is transient rather than steady. Compression 
and heating, i.e. the slow shock, is thus a primary effect of releasing magnetic energy, while 
the fast magnetosoic termination shock is a secondary effect caused by abruptly halting the 
energy release. We propose here that LTSs are manifestations of the primary effect: slow 
shocks. 

We present our model by first reviewing (in §2) how magnetic energy is converted to 
heat when fast reconnection occurs at a current sheet. We show that similar predictions are 
made by two-dimensional steady-state models, such as Petschek's, as by transient patchy 
reconnection models. The latter can be modeled as a simple shock tube problem. In §3 
we use the shock-tube problem to demonstrate how thermal conduction leads to density 
enhancement far beyond what Rankine-Hugoniot relations predict. We find a simple relation 
between the shock-tube mach number and the maximum density enhancement. We next use 
thin flux tube models to investigate how much density enhancement could occur within a 
reconnecting current sheet. In §5 we apply these single-fiux tube results to a flare consisting 
of numerous flux transfer episodes. We show how the EM of the SH-LTS scales with the 
rate of flux transfer. The constant of proportionality is determined by the angle between 
field lines reconnected at the current sheet. Section 6 then considers how the energy released 
by reconnection drives chromospheric evaporation to generate a second, cooler component 
of the plasma. Thus while the chromosphere contributes most the the EM in a flare, it is 
not responsible for the LTS. 



2. Density increase from reconnection shocks 
2.1. The equilibrium current sheet 



Fast mag netic reconnection in models such as that of iPetschekl (1196 



investigators (JForbes fc Priestlll983l : lYokoyama fc Shibatalll997l : iBirn et al. 



, or subsequent 



20011 ). occurs at 



a discontinuity in the magnetic field called a current sheet (CS). Figure [T ] shows a simple 
example of such a CS in a two-dimensional quadrupolar magnetic field ( iPriest fc Raadu 
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19751 ) anchored to four photospheric source regions labeled A-D. The coronal field contains 
free energy because its field lines interconnect the sources differently than would the lowest- 
energy field, a potential field. Here A and B are connected by fiux in excess of that required 
by the potential field; the excess flux, A^, is shaded grey in the figure. The excess magnetic 
energy can be decreased by a reconnection electric field within the CS breaking an A-B (red 
dashed) and C-D (blue solid) field line to create new field lines connecting A to D and C 
to B. This will decrease the flux excess A\E' thereby decreasing the free energy. 




Fig. 1. — A current sheet whose reconnection releases energy. A quadrupolar field is anchored 
to sources. A, B, C and D. In the end-on view (left), an equilibrium current sheet extending 
between Y-points at P and Q, separates downward field lines (red dashed) connecting A to 
B from upward (blue solid) connecting C to D. Viewed from face-on (inset) the current 
sheet separates field lines differing by angle A6. The strength of the vertical [B^, blue) and 
horizontal {B^, green) field components along the vertical F-F' are plotted on the right. 
The magnitude (|B|, magenta) includes the guide field (By) component. 



When the two-dimensional equilibrium includes a field component in the ignorable di- 



rection, a guide field [By, magenta), the current sheet separates field fines which are not 
exactly anti-parallel but make a finite angle A6 (see "face-on" inset). Reconnection at point 
R creates a new field line AP'RP"D, the A-D line, and CQ'RQ"B as the C-B line. These 
are not the shortest possible lines between their end-points. For example, the section Q'RQ" 
(thick black) can be shortened to a straight segment Q'Q" , along the lower edge of the CS, 
by pulling the bend through the CS. In general such reconfiguration would be hindered by 
other field lines within the equilibrium. A current sheet, however, represents a "crack" in the 
magnetic field through which field lines may pass unimpeded. This unimpeded shortening 
is the b asis of all fast m agnetic reconnection scenarios, but was posed in an equivalent form 



first bv IPetschekl (119641 ). 



All other things being equal, shortening field lines will reduce the net magnetic energy 
given by the integral 



Wm = -J |B(x)|^d^x = -Jd^J\B\di= Jd^^ , (2) 

z>0 

where the outer integral in the final two expressions is over all field lines composing the 
corona, and the inner integral, 6Wm/S^, is over the individual field lines. To reduce the 
total magnetic energy one must decrease the contributions of all fiux tubes on average. This 
is done by either decreasing their lengths or their average field strengths, or both. Changing 
a field line from Q"RQ' to the straighter Q"Q' will clearly do the former. Since the field 
strength is minimum along the CS edge (Q) moving the field line there will also decrease 
the average field strength. Current sheet reconnection therefore liberates magnetic energy 
in both ways. 

Prior to any reconnection the current sheet is part of a stable magnetostatic equilibrium 
(stable to ideal perturbations, which do not change magnetic field line topology). This means 
that its magnetic energy is a minimum and thus its field lines are, in a collective sense, as 
short as possible: AP'Q"B represents the shortest possible path between A and B for that 
particular field line. Shortening that field line would require deforming and lengthening 
so many neighboring field lines that the net energy would increase. It is clear that the 
central segments, P'Q" and Q'P" ., are relatively straight paths between the end-points when 
compared to the bends, P'RP" and Q'RQ", resulting from reconnection. 

That shortening a magnetic field line can also decrease its field strength is not obvi- 
ous and warrants further comment. The coronal magnetic field strength tends to decrease 
with height, so one expects retracting downward would increase rather than decrease field 
strength. This is in fact the case for the arcade field below the CS {z < zq in Fig. [1]) and 
probably describes the later phase of a fiare. It is therefore unclear whether such loops were 
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themselves releasing magnetic energy (their length decrease overwhelming the field increase) 
or were consuming energy being released elsewhere. 

The situation is very different within the CS, zq < z < zp. The tips of the current 
sheet, Q and P, are Y-type null points where the field components perpendicular the the 
current (i.e. the reconnection components shown in green and blue on the right of Fig. [1]) 
vanish. The magnetic field strength is thus a minimum at each tip: |B| = Bg. Within the 
region, zq < z < z^^, the field strength increases with height making it possible to both 
retract and weaken the field simultaneously. In addition to its more favorable field-strength 
profile, the region within the CS offers the most favorable fractional decrease in length: from 
Q'RQ" to Q'Q"; there appears to be less potential for shortening the legs of the post-flare 
loops, segments CQ' and Q"B. These factors combine to make the CS the region where 
most of the magnetic energy conversion will occur. All viable reconnection models are set 
within the CS and this is the region to which we will conflne further consideration. 

Shortening a flux tube might also decrease its volume thereby providing a means of 
raising its average density. The volume of a flux tube, per unit flux, is 



6V_ 



IBI 



(3) 



Decreases in length and fleld strength contribute in opposite senses and may therefore lead 
to either an increase or decrease in total volume as a result of energy release. Moreover, a 
flux tube anchored in the chromosphere is coupled to a reservoir from which it may change 
its own mass independent of the volume change. This is how evaporation can increase the 
density in the flux tube over long time scales. 

The main causes for density increase, at least initially, are slow magnetosonic shocks 
which are dynamical consequences of very rapid fleld line retraction. Dynamical evolution 
of shortening flux tubes has been treated using several seer aingly different models. In the 
earliest of these, the steady-state, two-dimensional model of iPetschekl (Il964f ). the reconnec- 
tion electric fleld occurs along the inflnite line z = zr wi t hin a current sheet without guide 
field {A9 = 180°). The model was generalized by ISowardI (Il982l ) to include a guide field and 
ther eby admit A^ < 180°. In a further generalization studied by iLinton fc Longcopd (120061 ) 
and lLongcope et al.l (120091 ). the electric field is localized in the erstwhile ignorable direction 
{y). A transient instance of such "patchy" reconnection produces two isolated fiux tube, such 
as the black and violet curves, Q'RQ" and P'RP" respectively, shown in Fig. [T] These were 
found to evolve approximately as thin fiux tubes, subject to their own magnetic tension. 
The layers of fiux separated by the CS confine the fiux tubes but otherwise let them evolve 
undisturbed. The density evolution turns out to be very similar in the two-dimensional, 
steady state model and the three-dimensional, transient model. 
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2.2. Thin flux tube retraction 



We illustrate this generic post-reconnection dynamics with solutions of thin flux tube 
equations within a CS. To simplify the geometry we ass ume a symmetric current sheet of 
the Green-Syrovatskii type (lGreenlll965l : ISyrovatskiilll97l[ ) between zq = and zp = h. The 
vertical fleld component (also called the reconnection component) varies across the sheet as 



B.{z) 



± 



2B 



zO 



h 



^Mh 



(4) 



where the upper and lower sign correspond to the front {x > 0) and back (x < 0) of the 
sheet respectively. There is a uniform guide fleld, By = Bg. We solve only for that portion 
of the fleld line within the CS (the Q'Q" segment) and only until effects of reconnection 
have reached these ends. We thereby restrict ourselves to times earlier than effects from the 
boundary, such as chromospheric evaporation, could play a role. 



The dynamical evolution is treated using the thin flux tub e dyna mics of lLinton &: Longcope 



( l2006l ). includi ng parallel dynamics added by iLongcope et al.l (120091 ) and fleld-aligned trans- 
port added by lGuidoni fc Longcopd (|2010[ ). The axis of the tube is a space curve. 



r(£,t) = y{i,t)y + z{i,t)z 



parameterized by length £. The unit tangent vector b = dr/di, is parallel to the internal mag- 
netic fleld. The low-/9 tube moves through the current sheet governed by its own dynamics. 
The velocity of an element changes due to internal gas pressure gradi ents, magnetic curva- 
ture, external magnetic pressure gradients and viscosity according to ( jGuidoni fc Longcope 
20111 ). but gravitational forces are neglected. 



P 



dt 



- dp B"^ dh 



477 



V±B' + B 



d 



di \B 



^bb 



dv 
'di 



(5) 



Here V^ is the component of the gradient perpendicular to b. While evolution of the curve 
changes the fleld's direction, its strength is set by pressure balance across the tube. The 
external pressure is provided by the flux outside the current sheet making it a flxed function 
of position 



B{y,z) 



Bl + Bl 



Bl 



'^'A 



< z < h 



(6) 



The tube is assumed to be small enough not to affect the sheet's equilibrium by its motion. 
The dynamic viscosity /i depends strongly on temperature T = {m/k^)p/p, where m is the 
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mean mass per particle and fee is Boltzmann's constant. Its temperature dependence is given 
by the coUisional Spitzer-Harm form 

/i ~ 0.12 dyne- sec (T/10^K)5/2 _ (7) 

This is small enough that it is often ignored, but it is essential for resolving the shocks which 
develop within the flux tube. 



We solve these equations using a Lagrangian numerical code (DEFT, iGuidoni fc Longcope 



2OIOI ) which follows the tube's mass elements and finds the mass density p using the elon- 



gation between pairs of elements and the cross-sectional area 5 A = S^/B. The resulting 



density undergoes an evolution satisfying ( JGuidoni fc Longcopell201ll ) 



dlnp dlnB 9(b ■ v) dh 

dt ^ ~di dr~ ^ ^"di ' ^^^ 

The first term on the right hand side (rhs) reflects the density enhancement from lateral 
compression as the tube moves into stronger magnetic field: p oc 5; it will move into 
stronger field if the reconnection site is above the sheet's midpoint {z^, > h/2). The second 
term reflects axial compression either from shortening or from shock compression. The third 
produces compression when radius of curvature is decreased through perpendicular motion. 

The pressure of a given mass element evolves according to an energy equation 
dp 5 dlnp 2 /^ dvV 2^d f k dT\ 

The first term on the rhs reflects adiabatic work, with 7 = 5/3, the second is viscous 
heating and the third is thermal conductiono The thermal conductivity, k, has temperature 
dependence identical to the viscosity and can be expressed as k = {k^/ffi)p/'Pi, where 
Pr = 0.012 is the Prandtl number using Spitzer collisional transporto Since we consider 
only a flux tube, all transport processes are automatically field aligned; we consider no 
viscous or thermal transport across the magnetic field. Radiative losses are neglected here 
because they are negligible at the high temperatures and short times scales of interest; we 
return to the point below. 



3 



Since heating occurs relatively slowly through viscosity the ther mal conduction is generally s mall and 



we find it unnecessary to use any form of conductive flux limiter (see iLongcope fc Bradshawll201Cll ) . 

^Thc Prantdl number is very small because viscous stress is transmitted through random ion motion, 
while most thermal energy is carried by the lighter, faster electrons (in this single-fluid treatment, a.k.a. 
MHD, both species have the same temperature) 
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The flux tube evolution begins after the reconnection electric field has created a bent flux 
tube, Q'RQ" in Fig. [2l (The other section, P'RP" evolves according to identical dynamics, 
so we disregard it.) The equations are initialized with uniform density and pressure and a 
configuration, r(£, 0), tracing curve Q'RQ". The reconnection electric field can in principle 
dissipate energy, heat ions or electrons, or even accelerate particles. All these possible effects 
are, however, ignored when the tube is initialized with uniform pressure and density. We 
choose to neglect these effects on the premise that the dissipation occurs within such a 
small volume that it cannot have a significant energetic effect on the subsequent dynamics. 
This same approximation, neglecting direct dissipation such as Ohmic heating, was used by 
Petschek and subsequent investigators of fast magnetic reconnection. 

Figure [2] shows the numerical solution of Eqs. (EJ, ([H]), and ([2]), by DEFT. The initial 
bend at point R decomposes into two rotational dis continuities, RD and R D' , and two gas- 



dynamic shocks, GDS and CDS', as described in iLongcope et al.l ( l2009l ). The rotational 
discontinuities propagate along the erstwhile equilibrium field lines at the local Alfven speed. 
They change the direction of the magnetic field, but not its strength or any plasma properties. 
Tension at the rotational discontinuity accelerates plasma along the bisector moving the 
plasma both downward {vz in orange, this is the retraction) and inward {vy in green). The 
sect ion between them is slightl y curved owing to the variation in this speed within the CS 



see 



Guidoni &: Longcopdl201ll ). but is basically horizontal. Thus the inward flow is actually 
parallel to the axis of the flux tube within the RD-RD' segment. While the Lorentz force 
(i.e. tension) accelerates only perpendicular to the local field, subsequent changes to the field 
direction renders the flow parallel. 

The inward parallel flow generated by the rotational discontinuities resembles a hydro- 
dynamic shock tube. The field bends to approximately horizontal through an angle A0/2. 
The field strength |B| is the same as in the equilibrium field and the plasma properties, in- 
cluding initial plasma beta, (3q, are those of the pre-reconnection ambient s tate. The parallel 
flow thus has a hydrodynamic Mach numbeio ( Longcope et al.ll2009l . l2010l ) 



/ 24 
M, = ^— sin2(A^/4) . (10) 

The solution in Fig. [2] begins with reconnection at A9r = 106° and (3or = 0.02, which would 
generate an inflow at Mi = 3.1. In the middle of the CS, z = h/2, where the field is stronger 
(/3o = 0.014) and the angle is greater {A9 = 120°) the flow could reach Mi = 4.6. By the 



""This should be confused with the Alfven Mach number often used in the reconnection Uterature to 
characterize the reconnection rate. M, is normahzed to the sound speed Cg and characterizes the speed of 
field line retraction following reconnection. 
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Fig. 2. — A thin flux tube simulation of post-reconnection relaxation. The initial field 
lines are shown in the upper left panel as red and blue dashed curves. These reconnect 
at point R (asterisk) and relax to the solid black curve at a later time. Fluid within the 
tube moves downward {vz orange) and inward {vy green). The profiles of these velocity 
components are shown as curves below, normalized to Cso, and the peak of Vy, corresponding 
to Mi, is represented as a green dashed horizontal line. Both velocity components are 
initiated by magnetic tension at the rotational discontinuity, RD' (square). The horizontal 
velocity is reduced at the gas dynamic shock, GDS' (diamond). The density (p, magenta) and 
temperature (T, blue) are plotted below this, normalized to their ambient values. Horizontal 
dashed lines of the same color show the Rankine-Hugoniot values corresponding to the Mi 
from above. The density changes at the slow shock; the temperature increases within the 
heat fronts (HF, triangles). The central value (xs) undergo a time-evolution plotted in the 
lower right panel against height, z (increasing leftward). 
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later time shown in the figure the actual parallel fiow has a Mach number Mi = 4.2 (see 
middle panel). 

These infiows should not be confused with the slow "infiows" feeding the reconnection in 
steady state reconnection models. In the more conventional end-on view, these appear fiowing 
toward the current sheet, where they are accelerated, by rotational discontinuities, to form 
the "outfiow jet", directed away from the reconnection site (upward and downward). The 
infiows shown in Fig. [2] are actually part of the same "outfiow jet", having been accelerated 
by the rotational discontinuities. Since they are part of the outfiow jet, they are also present 
in steady-state models with a guide field, but the component directed along the ignorable 
direction (i.e. the infiow, Vy) is not often discussed. In the less- conventional face-on view, 
shown in the figure it is evident that these fiows are natural consequences of field line 
shortening. 



2.3. Shocks 

The supersonic, inward, parallel fiows from each rotational discontinuity collide in the 
center creating gas dynamics shocks where the parallel fiow speed drops to zero (green curve 
in middle panel of Fig. [2]). Both the density and temperature increase across such a shock, 
resulting in a hot dense central plug (magenta and blue curves in the bottom panel). 

The gas dynamic shocks differ in several significant respects from the slow magne- 
tosonic shocks of the classical Petschek model, especially its anti-parallel version {A9 = 180°, 



Longcope et al.l 120091 ). When the two-dimensional steady model includes a guide field, how- 



ever, there are two sets of shocks: rotational discontinuities enclosing an outfiow jet, inside of 



which are slow mag netosonic shocks where heating occurs ( lSowardlll982l : iForbes et al.lll989 



Skender et al.ll2003l ). For cases of non-trivial guide field (A^ < 150°) and low /3, the slow 
shocks are very nearly parallel shocks resembling gas dynamic shocks. Figure |3] shows only 
a modest decrease in magnetic field strength {B2/B1 ^ 1) and very small change in field 
direction, fii 2, across the slow shock. These cases look far less like Petschek's original model 
with a hot outfiow jet, than like the thin fiux tube model shown in Fig. [2l The greatest 
difference is a superficial one: a steady outfiow jet (Petschek) is replaced by a single plug 
moving downward and expanding outward. The latter becomes more akin to the former 
when multiple reconnection events are combined to generate a sequence of plugs. 

Many applications of reconnection models, either Petschek's or the patchy thin fiux 
tube model, predict the density and temperature enhancement from shock jump condi- 
tions, sometimes called Rankine-Hugoniot relations. Such relations follow from conservation 
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Fig. 3. — Jump conditions across steady shocks for different reconnection angles A6', accord- 
ing to different models of post-reconnection dynamics. The top panel shows the plasma /3 of 
the post-shock material, and the r niddle shows the density ratio. Solid curves are from the 
2d steady state m odel of ISowardl ( 119821 ). the diamonds are for the Id Riemann prob lem of 
Lin fc Led (Il994j ) and the dashed for the thin flux tube model of iLongcope et al.l (|2009[ ). The 
latter corresponds to a shock tube for inflow Mach number given by Eq. (TTOjl . plotted along 
the top axis. A x marks the conditions at point R of the run shown in Fig. |2l The bottom 
panel shows the jump in magnetic field strength (against the left axis) and the change in 
field direction (against the right axis) for the two MHD models. In the thin flux tube these 
are fixed to be B2/B1 = 1 and f2i_2 = 0. 



l aws and from the assumpti on that the shock is steady when viewed in a co-moving frame 
( ICourant fc Friedrichslll948l ). Under these conditions the density, temperature and magnetic 
field strength jump from an upstream value (subscript 1) to a downstream value (subscript 
2). The ratio of downstream to upstream values depends on the Mach number Mj which 
follows in turn from the angle A6. This angle is a parameter of the pre-reconnection current 
sheet in both thin flux tube models and steady, two-dimensional models; it should not be 
confused with the opening angle of the outflow jets when the latter solutions are viewed 
edge-on. 



Figure |3] shows the density ratio predicted by Rankine-Hugoniot relations for two- 
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dimensional steady state reconnection (solid) and thin flux tube dynamics (dashed). A 
third version of the r elations, obtained from a one-dimensional Riemann problem beginning 
with bent field lines (JLin &: Leelll994j ). is plotted using diamonds. This simplified problem 
captures the non-resistive evolution following reconnectio n and has been s hown to reproduce 
properties of two-dimensional steady state reconnection (JLin fc Led 119991 ). 



8M2 + 4Mi^4M2 + 9 + 9 



;ii) 



The Rankine-Hugoniot density ratio for the shock tube (dashed curve in the middle 
pane of Fig. [3]), 

Prh ^ P2 ^ 

Po Pi 2Mf + Mi^/IM[T^ + 9 

has a well-known asymptote P2/P1 — )■ 4 at large Mi (i.e. A6' — )• 180° and /3o — )■ 0). In 
this same limit the slow shocks of the MHD models (2d steady and Id Reimann problem) 
approach switch-off shocks — the limit originally used by Petschek. A switch-off shock 
annihilates the ma gnetic field {B2/B1 — )■ 0) and has a limiting density ratio P2/P1 -^ 2.5 
(JForbes et al.lll989l ). For less severe angles, say A6 < 120°, the slow shock barely defiects the 
magnetic field, (^1,2 < 3°) and therefore resembles a gas dynamics shock. Since this is the 
same as in the shock tube all three models show very similar behavior in this regime. There 
is still some decrease in field strength across the slow shock which is naturally absent from 
the shock tube. The former thus has a slightly lower density enhancement than the latter, 
although both are above the switch-off limit of 2.5. The very large post-shock /3 (top panel) 
is dominated by plasma pressure increase since the magnetic pressure decrease is small (it is 
zero for the shock tube). 

The density enhancement in the time-dependent solution does not, however, match 
the steady Rankine-Hugoniot value of Eq. (fTTj) . The lower panel of Fig. [2] contains dashed 
horizontal lines showing the jumps predicted by Rankine-Hugoniot conditions for Mj = 4.2 
(the observed Mach number). The central density notably exceeds this value; in fact it 
exceeds the maximum possible ratio of 4. This is due to a "transient" density overshoot, and 
temperature undershoot, which has been obser ved in previous time-dependent thin flux tube 
simulations (JGuidoni fc Longcopell2010l 1201 ll ). During the earliest phase of the retraction 
(see lower right panel of Fig. |2]) there is a brief period where the density ratio reaches 10: 
far above the theoretical maximum of Rankine-Hugoniot shocks. 

This density overshoot is a direct consequence of transient behavior resulting from ther- 
mal conduction. When transport is dominated by thermal conduction (i.e. when Pr <^ 1) 
the s t eady shock proflle include s a heat front ahead of the primary density jump ( iThomas 



1944J: lGradlll95ll: lGerma in"l96o[ lKennell[l988[ IPorbes fc MalherbelEgsel : IXu fc Forbes 



1992 



Yokoyama fc Shibatalll997i) . These are evident in the temperature proflles of Fig. [21 where 
the leading edge is marked by a triangle. For typical coronal parameters the steady state size 
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of the heat front is so large that is cannot be achieved during retraction ( IGuidoni fc Longcope 
20101). Thus the entire relaxation process occurs out of steady state as the heat front con- 
tinues expanding toward its asymptotic size. 

It seems that the actual density enhancement caused by reconnection shocks follows from 
transient shock evolution in the presence of thermal conduction rather than from the steady- 
state Rankine-Hugoniot values. The actual enhancement can be several times larger than 
the steady Rankine-Hugoniot values. This is a significant departure from standard practice 
where Rankine-Hugoniot iur np conditions have been widely used to estimate emission frorn 



fast magnetic reconn ection (jPetschekl Il964l : ICargill &: PriestI Il982l : IVrsnak fc Skenderl 12005 



Longcope et al.ll2010l ). While we have observed the discrepancy in thin flux tube models, 
these closely match MHD models in other respects so we hope to apply thin flux tube results 
generally. 



Thermal conduction transients in shocks 



The density overshoot can be studied thoroughly in a simplified, one-dimensional shock 
tube problem. Equations of continuity, momentum and energy conservation, along the linear 
tube, parameterized by length s, are 
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(14) 



using 7 = 5/3 as before. Viscosity p and conductivity k are computed using Spitzer-Harm 
form, Eq. ([7]), as they are in the thin flux tube model. This simple set of equations is 
equivalent to Eqs. ([5]), (IHl) and (l9l), for a static, straight flux tube (£ — )■ s and b ■ v — )■ m). 

The equations are solved within the interval, < s < L, with a rigid wall at the left 
boundary (s = 0): u = and dp/ds = dp/ds = 0. The wall represents the central stagnation 
point in the horizontal segment [RD-RD') of the thin flux tube solutions, and values there 
will henceforth be referred to as "central" . A steady, uniform inflow is introduced at the right 
boundary, u{L) = —Mi Cso, p(0) = po and p(0) = po? where Cso = yTPo/po is the ambient 
sound speed and Mj is the inflow Mach number. This represents the rotational discontinuity 
at which the parallel flow is created and its location is irrelevant for the gas dynamic shock; 
in practice it is repositioned further right whenever an effect from the wall approaches it. 
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The natural length scale in the problem 

«.„ = Jf^ = 43km (-§-V (j^p^y . (15) 

is related to the ion mean free path. The corresponding time scale, related to the ion collision 
time, is 

Tin = — = 0.25 sec — - — — . 16 

Cso \WK) VlO^cm-^/ ^ ' 

The system is initialized with uniform properties, p(s, 0) = po and p(s, 0) = po- The 
velocity is initialized with a smooth step, 

m(s,0) = -Mitanh(s/A) , (17) 

whose gradient scale is denoted A. It is necessary to position the right boundary, L S> A, in 
order that the initial condition be consistent with the boundary condition. 

The solutions to these equations, such as the example shown in Fig. HJ show features 
similar to the retracting flux tube of Fig. |2l Density piles up against the wall (s = 0) until 
a pressure is achieved capable of driving a shock into the flow. This compression raises the 
temperature of the high-density plug, but thermal conduction leads to extended heat fronts 
(triangles) in front of the shock. The pressure achieves its flnal (Rankine-Hugoniot) value 
promptly, and thereafter remains fairly constant. The peak temperature rises more slowly 
due to thermal diffusion spreading the heat into the broad region behind the front. The 
density therefore overshoots its flnal value and decreases thereafter. The fully developed 
conduction front is extremely large, and is therefore achieved only after a long time. As it 
is approached the temperature climbs and density falls - both rather slowly. 

Since the tube length L is irrelevant, /i deflnes the length scale £j05 and the Prantdl 
number is flxed (Pr = 0.012) the only two genuinely free parameters of this system are the 
Mach number, Mj, and the initial shear length A. Together they deflne the ballistic collapse 
time, Tb = Tjo^io/AMj, taken for the initial velocity shear to steepen into a shock. 

The size of the initial velocity gradient. A, relative to £jO) determines two slightly different 
patterns in the evolution of the central density (see Fig. E]). For small initial gradients, 
A ^ ^jO) the viscous heat generated during the steepening is more than compensated by the 
thermal conduction. As a consequence the central entropy decreases slightly during shock 
formation, causing p(0) to evolve sub-adiabatically. In the opposite limit, A <^ ^jo, the initial 
velocity shear is rapidly dissipated by viscosity before the density can respond. This leads to 
a rapid increase in central entropy which is later offset by thermal conduction. The pressure 
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Fig. 4. — The evolution of the shock tube solution for Mi = 2.0, A = lOiiQ. The right column 
shows the profiles, p{s,t) (top), p{s,t) (middle), and T{s,t) (bottom) at a sequence of times. 
Successive profiles are displaced rightward, and symbols show s = for each profile. An axis 
along the bottom of each shows s/iio at t = 0, while a second along the top of the upper 
plot shows the same coordinates as they pertain to the final (right-most) time (t = IOOtjo). 
As on the bottom plot show the extent of the heat front. The left column shows the central 
value, p(0,t), p{0,t), and T(0,t), over time (scaled to Tjo). Symbols correspond the times of 
the profiles on the right. The vertical dotted line shows the time of peak density. Horizontal 
dashed lines on all panels show the Rankine-Hugoniot values for each quantity. 
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Fig. 5. — Density evolution for cases of Mj = 2 (bottom) and Mj = 5 (top). Left panels 
show the central density, p(0,t), versus time for initial conditions characterized by different 
values of A/£jo- Triangles denote point of peak central density. In the lower left panel these 
values are (reading peak values clockwise from the lowest) X/to) = 0.1, 0.3, 1, 2.5, 5, 7.5, 
10, 15 and 25; in the upper left panel they are \/iio = 25, 40, 60, 100 and 150. The right 
panels show the same evolution plotted as p(0,t) vs. T(0,t). Vertical and horizontal dashed 
curves are the Rankine-Hugoniot values. The dotted curve dividing grey from white regions 
is the adiabatic curve through the initial state; the grey region has entropy below initial. 
The other dotted curve is the locus with pressure matching the Rankine-Hugoniot state. 
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from this initial heating hinders the later collapse, so cases with A ^ fjo, tend to reach lower 
peak central densities than do the smaller initial gradients. 

The density overshoot is a consequence of the initial growth of the heat front by thermal 
conduction. Except for cases of steep initial gradients (A <^ £jo) the initial evolution is nearly 
adiabatic, and then becomes slightly sub-adiabatic as conduction cools the central fluid, as 
just described. The collapse time, Xb, increases with A causing peak density to occur later, 
as shown in Fig. O The peak value achieved is, however, little different than for modest 
values of A. For low mach number of Mj = 2 (bottom row) the largest density enhancement 
occurs for A = lOiio, at t = 9 ^ 2rb. For the large Mach number, Mj = 5 (top row), the 
peak is at A = lOOiio, at t = 27 ~ 1.5rb. 

The initial conduction is so effective at cooling the central plasma that the peak density 
occurs just above the intersection of the adiabatic curve, p ~ p^^^, and the final (Rankine- 
Hugoniot) pressure p^-^^. This point is at the crossing of the two dotted curves on the right 
panels of Fig. [51 The adiabatic density enhancement 

, 3/5 
Pad (PnH^ 



Po V Po 

is thus a lower bound to the maximum density overshoot. Figure [6] shows how this en- 
hancement increases rapidly with Mj. For Mj > 1 this is significantly greater than the 
steady-state (Rankine-Hugoniot) enhancement (dashed) which asymptotes to 4. The maxi- 
mum enhancement for runs with optimal A ( x ) lie above the adiabatic curve. An empirical 
fit Pmx = 1-33 Pad (dotted) fit the data well. 



4. Post-reconnection enhancements 

The thin flux tube model of post-reconnection retraction differs from the shock tube 
problem in several respects. First it is initialized at rest, rather than with a prescribed 
velocity profile like Eq. flTT]) . The inflow is generated not by a right hand boundary, but at 
a rotational discontinuity moving horizontally away from the center at some fraction of the 
local Alfven speed. The rotational discontinuity travels along a field line with varying field 
strength and direction. This variation leads, in turn, to variation in /3o and A^, and thus 
in Mj according to Eq. (TTOj) . The gas dynamic shocks will therefore come from inflow with 
time-dependent properties. Finally, the variation in field strength causes its own variation 
in density acco rding to the first terrn on the rhs of Eq. ([H]); this latter effect was extensively 



investigated by lGuidoni &: Longcopd (l201ll ) 



In spite of these differences, the central density enhancement and heat fronts of the 
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Fig. 6. — Peak density enhancements, niax[p(0, t)/po]; for shock tube solutions with mach 
numbers Mj. Crosses are solutions with optimal gradient scale A. The solid curve shows the 
adiabatic enhancement from Eq. ( ITSl) . Dashed curves are the Rankine-Hugoniot enhance- 
ments from Eq. fITT]) . Left panel shows the enhancement vs. Mach number; right panel shows 
it vs. the temperature at time of peak density. 



thin flux tube simulations are very similar to those of the shock tube. This is evident by 
comparing profiles on the right of Fig. S] to the bottom panel of Fig. |2J Both cases have 
heat fronts far ahead of the density jump (slow shock). In both cases the central density 
overshoots its Rankine-Hugoniot value and then approaches it from above, as a direct result 
of the growing heat fronts. 

Figure [7] shows several other examples of central density enhancements from thin flux 
tube solutions. All occur within identical current sheets: A^q = 120°, and height h = 
3 X 10^ iio- In each case the reconnection is initiated at different positions zji and the Pr was 



adjusted to keep the initial value of M, = 7, according to Eq. (jTOl) . In all solutions the central 
density rises rapidly to a value above the Rankine-Hugoniot limit (the horizontal dashed line 
on the lower left). In reconnections closer to the bottom the density falls immediately after 
its initial rise. For the highest reconnection point, zn = 0.875/i (diamond), the progression 
into stronger fields and greater angles keeps the density ratio at p/po — 13 — far above the 
Rankine-Hugoniot limit. 
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Fig. 7. — Density evolution for thin flux tube solutions in identical current sheets with 
reconnection initiated at positions, zr/H = 0.875 (diamond), 0.75 (x), 0.5 (triangle), and 
0.25 (square). The field strength at the reconnection point is plotted on the upper right curve. 
Value oi /3ji was takes tohe /3r = 6x 10~^ (diamond), 5 x 10"^ (x and square) and (3 x 10"^) 
in order that the initial Mach number was Mj = 7. The bottom row shows the evolution of 
the peak density versus time (left) and height (right; note that z decreases rightward, which 
is the sense of the evolution). The upper left column shows the total emission measure of 
heated plasma versus time. The dashed curve shows the result for a steady Mj = 7 shock. 
The bottom time axis is scaled to Xjo; the top axis uses Ugo = 3 x lO^cm"^ and Tq = 2 MK 
to convert this to seconds. 
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When plotted agains t/xjo (lower left) the central density shows its further resemblance 
to the shock tube. Since most reach their peak before t = IOtjo, even though Mi = 7, these 
would seem to resemble shock tube solutions with A <ti iio- This would explain their tendency 
to reach a peak value below the adiabatic, p = Po(Prh/Po)°'^) shown as a horizontal dashed 
curve. Other runs reveal that reducing £io, and thus Tjo, has an effect similar to increasing A: 
peak density is raised still further. Doing so at high Mach number is numerically challenging 
due to the very small scale created at the center. 

Additional comparison with observations can be made using the emission measure EM 
of the reconnection-heated plasma. Thin flux tube solutions yield EM per unit flux in 
analogy to the field line integral in Eq. ([2]), 

6EM r ^di 

nl-B , (19) 



5$ J ^ B 

where this integral is only over the heated portion of the tube (i.e. where T > To). Time evo- 
lution of this quantity is plotted in the upper left panel of Fig. [3, normahzed to nl^h/Bg. In 
each case the EM continues to increase after the density has peaked, owing to the increasing 
mass swept up by the shock and heat front. 

The EM begins dropping before the simulation ends, as the heat front and rotational 
discontinuities generate a g rowing low-density void (rarefaction wave) ahead of the shock 



( JGuidoni fc Longcopell201ll ). The dashed line in the upper left of Fig. [7| sho ws the emission 



measu re predicted by the the simple Rankine-Hugoniot analysis used by iLongcope et al. 



( 120101 ): density and temperature jumping at a steady jump moving at the Rankine-Hugoniot 
shock speed. All cases exceed this prediction due to their much higher density as well as the 
greater extent of their heat fronts. 



4.1. General behavior 

The general behavior of the thin flux tube model can be gleaned from a composite 
of many different solutions. The bottom panel of Fig. [8] shows the maximum density en- 
hancment vs. Mi for 25 different runs. The collection spans central reconnection angles 
90° < A^o ^ 133°, mean free paths 7 x 10^^ < iio/h < 7 x 10^'^, and reconnection points 
0.75 < z^ < 0.82 except for the examples shown in Fig. [7] (indicated by symbols). Runs 
with the highest Mach numbers terminate before complete retraction due to the arrival of the 
heat fronts at the loop footpoints (this is discussed further below). This is why plots against 
z in Fig. [7] (lower right) end short of the origin. A more systematic sample of runs would 
clearly be useful, and will be attempted in a future investigation. The present hodge-podge 
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does, however, serve our present purpose by illustrating how typical solutions compare to 
the Rankine-Hugoniot model, plotted as a dashed line. 
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Fig. 8. — Thin flux tube solutions over a range of Mach numbers, Mj. The 4 runs from 
Fig. [7] are designated by the same symbols (diamond, x, triangle and square); the run from 
the following section is designated by a *; all others are plotted by a +. The top, middle 
and bottom panels plot 5EM/5^, from Eq. (1201) . temperature and density enhancement 
respectively. The dashed lines correspond to the Rankine-Hugoniot limit from Eqs. fl22]) and 
( ITTj) in the top and bottom panels. The broken line in the bottom panel is the adiabatic 
density ratio given by expression ( !T8|) . The symbols on the top panel given 6EM/6^, when 
the simulation ends, often when the heat fronts reach Q' or Q". Vertical upward lines estimate 
how high it would be if the simulation could continue until the retraction was complete. 

The density enhancements in Fig. |8] (bottom panel) lie well above the Rankine-Hugoniot 
value (dashed) in all the myriad runs. The latter model of density enhancement appears to 
be needlessly conservative. The cases with higher Mach numbers tend to have enhancements 
below the adiabatic limit (broken line, Eq. [IS])- This is probably due, as remarked above, to 
the similarity of the thin flux tube solutions to shock tubes solutions with A ^ £jo- Density 
ratios at or above a factor of 10 appear to be the rule, with exceptions most often for low 
reconnection points (square) or extremely small CS (i.e. large iio/h). 



The maximum observed temperature ratio, plotted in the middle panel of Fig. |S1 is 
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approximately limited by the Rankine-Hugoniot value (dashed curve). The runs at highest 
Mach number require the longest to reach steady shock conditions (i.e. Rankine-Hugoniot 
assumptions), and are the ones stopped before complete retraction. As a result their tem- 
perature maxima tend to fall well below the Rankine-Hugoniot value. Even in these cases, 
however, the peak temperature is 10-20 times the starting value, easily eligible for a SH 
classification. 



4.2. Mean properties: super-position 

The results above apply to a single episode of one post-reconnection flux tube retracting. 
We assume a solar flare consists of hundreds of such episodes as flux is transferred through 
the current sheet at hundreds of separate patches. The net properties of the flare follow from 
a super-position of episodes resembling those in Fig. [71 If the flux transfer occurs at a mean 
rate $, the mean EM of the super- hot source will be 

,^_, ■ f SEM , . 5EM , , 

where the time integral is over a single episode. Time integrals of the four curves in Fig. 
[7] yield bEMjb^ = 0.03, 0.09, 0.14 and 0.10 times the dimensional factor nlf^h^/csoBg, for 
zn/h = 0.25, 0.5, 0.75, 0.875. Even though the highest reconnection point {zr = 0.875 /i) 
reaches the greatest density, its heat fronts reach the footpoints earliest giving it a shorter 
duration and thus a smaller net emission measure, SEM/6^, than the reconnection at zr = 
0.75 h. The dimensionalizing factor 

^^ = 6x10^^ ''' ( "-° Y( ^ Vf ^° V^^Y^. V^ .21) 

CsoBg cm3Mx Vl09cm-3; l^lOOMm/ VlO^Ky \10GJ ' ^ ' 

gives a reasonabl e LTS, (EM^^) '^ 10 ^^ cm"^, for reconnection at a typical mean rate $ = 
2 X lO^'^Mx/sec JLongcope et al.lboiol ). 



The top panel of Fig. [8] show 6EM/6^ for all 25 runs in the collection. In cases where 
the run terminates before full retraction, an attempt is made to correct the value to a full 
retraction. The instantaneous EM at the run's end, which is typically lower than the peak, 
is extended to the time the tube would have reached z = 0. This extrapolation is indicated 
by a vertical line upward from the value actually observed (symbol). Since faster heat fronts 
occur for larger Mj (the right of the plot), those are the cases with vertical extensions. 



The net EM can be compared to the prediction of a steady-shock model ( JLongcope et al. 
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2010|) 



'-^-&){t^M -'^^«^^' --'^^'^"^ > ■ '^^' 



plotted as a dashed curve on the top panel of Fig. [S]for zr = h/2 an d A6'b = 120°. This 



would be the EM in a steady 2.5-diniensional model of the iPetschekl ( 1l964l ) variety. The 
inverse scaling with Mi results from the faster retraction speed, relative to the normalizing 
speed, Cso, when f3ji becomes small (and thus Mj becomes large). The thin flux tube solutions 
have somewhat greater EM than the Rankine-Hugoniot model would predict, but appear 
to share its decrease with increasing M,. 



5. Energetics 

To illustrate the energetics of creating of the SH-LTS, and its subsequent evolution, we 
consider a specific CS suitable for a compact flare. The CS is h = 50 Mm tall with A^q = 120° 
and a guide field -B^ = 86 G (see Fig. [9]). Reconnection is initiated at zr = 0.75 h = 37.5 Mm 
where A6r = 113° and Br = 156 G. Properties of this run are plotted with an asterisk on 
Fig. [HI The right panels of [9] shows the heat fronts preceding the rotational discontinuities, 
and the upper left shows that they reach points Q' and Q" at t = 5 sec, well before the 
rotational discontinuities do. 

The initial tube, Q'RQ" is L{0) = 98 Mm long, and retraction shortens it to Lmin = 60 
Mm, the distance between Q' and Q". Remarkably, while the length is decreased by 40%, 
the volume of the flux tube remains virtually unchanged due to weakening of the field. Thus 
no density enhancement can be attributed to compression through field line retraction. 

The magnetic energy, per unit flux, inside the tube is given by the lin e integral in Eg. 



2]). The confining external field does work on the flux tube as it retracts (jLongcope et al. 



20091 ) ■ contributing additional energy to an effective potential 



^^ = — \B\de . (23) 

5$ 47ry ' ' ^ ' 

Since thin flux tube equations, Eqs. ([5])-([9]) are approximately conservative (to ~ /3o), any 
decrease in this potential energy appears as either thermal energy or bulk kinetic energy. 
In its initial state the flux tube has a potential energy SW^^ /5$ = 12 x 10^'^erg/Mx. 
Shortening and weakening reduces this by 66% to 

Eo = — rr^ = —^, = 4 X 10 erg/Mx . (24) 

o<P 47r 
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Fig. 9.— A run with Mj = 6 reconnected at zr = 0.75 /?, in a CS with A6'o = 120°. The 
central panel shows, for illustration purposes, an end-on view of the CS, with an asterisk 
marking the reconnection site. The face-on-view of the FT is plotted to the right of this; 
the dashed curve is the initial state, and sohd curve is a later time (t = 3.5 sec). Squares 
and triangles mark the heat fronts and rotational discontinuities respectively, and a diamond 
marks the tube's center. The temperature profile at this same time is plotted below it on 
the same scale. The panel to the left of the end-on perspective shows the vertical positions 
of the heat fronts (dashed with square) rotational discontinuities (broken with triangle) and 
tube center (solid with diamond) vs. time (in seconds). The symbols occur at the time of 
the solid curve in the upper right. Below this the various energies are plotted against the 
same time axis, normalized to the energy minimum, Sq, given in Eq. (^^. The solid curves 
are, W^^j , K±, Wt and i^y, from top to bottom. The broken line is Wt + K^^, which remains 
confined to the flux tube after retraction has terminated. 
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Retraction will therefore liberate 



A^ 



M 



5$ 



^0 = 8x 10'"erg/Mx (25) 



t=o 
energy available to the power the flare. 



The ratio of liberated energy to Eq is not equivalent to the more common ratio of free 
energy to the energy of the potential field. The equilibrium current sheet carries a net current 

h 

^cs _ ^ I T3 M - '^ ^-0 f rjT :y , B,oh _ „ ^ , „io 



— iB-dl = -^ [ Jzih - z) dz = ^^ = QxlO^OMx/cm , (26) 
An J n h J 8 



after using Eq. (^ for the vertical field at the CS. (Translating to SI gives Ics = 9 x 10^^ 
Amps, typical of pre-fiare currents). Denoting by A\E' the flux discrepancy giving rise to 
the current sheet, the flux which must be reconnected (see Fig. [1]), the free energy in the 
equilibrium can be derived from integrating the work required to establish the current sheet: 
AW = (J/c)A\Ef/2. If the present case has a typical pre-fiare flux, A\l/ = 3 x 10^^ Mx, its 
free energy would be AW = 3 x 10^^ ergs. The potential field of the entire active region will, 
however, include more flux, $^p ~ 3 x 10^^ Mx, and thus contain energy about an order of 
magnitude larger. Multiplying Eq. f l25|) by A^ gives twice AW because AS quantifies the 
differential energy release at the full CS size h. As h decreases so will AS; integrating the 
release over the entire reconnection, /i — )■ 0, approximately matches AW. 

To dimensionalize the plasma response we choose ambient density and temperature 
neo = 8 X 10^ cm-3 and Tq = 3 x 10'' K so that /3r = 0.007 and Mi = 6.0 according to 
Eq. ( [T0|) . The temperature profile at an early time (see the lower right of Fig. |9]) shows the 
heat fronts (squares) running far ahead of the rotational discontinuities (triangles). After 
5 seconds the fronts have reached the bottom of the CS (points Q' and Q") even though 
the center of the flux tube is only at the midpoint. To study evolution beyond this time 
we include additi onal straight segments run ning along the CS edge, for the heat fronts to 



propagate along ( IGuidoni fc Longcopell201ll ). This permits the solution to be followed until 
the center of the retracting flux tube (diamond) reaches the bottom of the CS (at t = 10 
sec). This occurs before the rotational discontinuities reach the bottom, so the tube is not 
quite straight: L = 1.07 Lmin- As a result the potential energy is still 22% above So (see Wm 
curve in Fig. [9]) . 

As the tube retracts it converts potential (magnetic) energy into thermal energy and 
kinetic (per unit flux), given by the field line integrals 

6Wt 3 f di 5K 1 /• , ,2^^ 

7^ = 2/^5 ' 5$ =2/^""" 5- ^27) 
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The kinetic energy can be further decomposed into that from motion perpendicular {K±) and 
parallel {K\\) to the tube's axis b. The lower left panel of Fig. [9] shows how all these energies 
increase as the potential energy decreases. The majority of the potential energy decrease, 
0.78A£, is converted into perpendicular motion: the rapid retraction. The remaining energy 
is converted initially into parallel flow {K\\) the collision of which generates the gas dynamic 
shocks which thermalize a fraction. By the end of this simulation 9% of the liberated energy 
has been thermalized while 13% remains in parallel flow. While the thermal energy is 
concentrated within the central mass concentration, the plasma (3 never exceeds 0.8. 



6. Density evolution following retraction 

The model above concerned only the initial stage in post-reconnection evolution: the 
reconnected flux tube retracting through the CS. To place this into the context of the entire 
flare we briefly discuss the subsequent stages illustrated in Fig. (TUl The downward retraction 
of the flux tube will be stopped by the closed fleld lines underlying the current sheet, the 
arcade (Fig. [TOb). At this point the fleld line has shortened as much as possible, and its 
magnetic energy has been converted mostly to bulk kinetic energy in Alfvenic downward 
motion, Kj_. This energy must somehow be dissipated if the flux tube is to join the static 
arcade. In a steady-state model this dissipation occurs at a standing fa st magnetosonic 



termination shock shown in the flgure (JForbes fc Priestlll983l : lForbeslll986l ). In a 3d, time- 



dependent model, such as the thin flux tube, it is more likely to be a traveling shock or a 
pulse of MHD waves initiated by the impulsive encounter between retracting flux tube and 
static arcade. After this occurs the static loop will enter a cooling phase. 

The different phases may be intermixed if heat fronts drive chromospheric evaporation 
during the retraction. Figure [9] shows the heat fronts reaching the edge of the CS {Q' and Q") 
far ahead of the tube itself. From there they would continue along the arcade legs of the tube 
eventually reaching the chromospheric footpoints. Conductive flux into the chromosphere 
would drive evaporative upflows into the loop. If the heat fronts were fast enough this upflow 
could reach the rotational discontinuities and modify the initial stage. This seems unlikely 
in the present case where the rotational discontinuities reach the CS bottom only 5 seconds 
after the heat front. We therefore neglect this possibility for the moment and return to 
discuss its implications below. We prevent any back-reaction from the tube footpoints by 
adding very long straight sections, as described above. 

Assuming the retraction occurs in the absence of chromospheric evaporation allows us 
to clearly distinguish that evolutionary phase from the subsequent cooling phase. Figs. [TOb- 
[TUb . In this later phase the loop has achieved magnetic equilibrium within the arcade. 
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Fig. 10. — Three stages in the post-reconnection relaxation within a current sheet viewed 
from the end-on perspective. Dark gray shows a particular flux tube at (a) t = just 
after reconnection at point z = zr-, (b) t = tb after being retracted through the current 
sheet, and (c) t = tc after additional flux (light grey) has been reconnected. The additional 
reconnection between tb and tc decreases the vertical extent, h, of the current sheet, moving 
point Q upward. The dark grey flux tube moves downward by Az (arrow) in a process known 
as "shrinkage". Its volume also decreases as its apex field strength increases: "collapsing 
trap" . 



The energy left in the loop from its retraction drives chromospheric evaporation by which 
radiative cooling is enhanced. It is this much longer phase that is observed in soft X-ray and 
EUV images as post-flare loops. Subsequent retraction episodes will pile additional flux on 
top at the same time it decreases the CS size h. The loop in question therefore "shrinks " 
downward, as several investigators have observed (JForbes fc Actonlll996l : iReeves et al.ll2008l ). 
This shrinking occurs relatively slowly, and the loop remains in magnetostatic equilibrium 
during the process. We therefo re neglect any energeti c effect it might have, including the 
"collapsing trap" contributions (ISomov fc Kosugil 119971 ). 



The energy powering the cooling phase comes primarily from the energy released during 
the retraction phase. The majority of this energy (almost 80% in this example) is in per- 
pendicular motion of the flux tube {Kj_). It is not, however, clear how much of this will be 
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dissipated within the flux tube since FMS waves propagate perpendicular to the field. With- 
out a clear means of resolving this point we restrict consideration to the remaining energy, 
0.22A£^, in parallel fiows and thermal energy. This fiux-tube-confined energy is plotted as a 
broken line in the lower left p anel of Fig. [HI Perpendicul ar forces from the underlying field 



will have little affect on this ( iGuidoni fc Longcopd 120111 ) so it will remain within the fiux 



tube even after it has come to rest. Eventually we expect all parallel fiow to be thermalized, 
leaving SWt/S^ = 1.8 x 10^°erg/Mx in the static tube. This is the energy powering the 
cooling phase. 

There are few truly simple models of chromospheric evaporation, especially driven by 
thermal conduction. Since our main interest lies in the evolution of temperature and den- 



sity within the loop, we use a zero-dim ensional mode . 



Antiochos fc Sturrockl ( 119781 ) and later by lCargill et al.l(ll995[ ). More recently iKlimchuk et al 



of t he kind origina lly proposed by 



( 20081 ) modified this treatment to account more accurately for the energetic contribution of 



chromospheric evaporation through the enthalpy fiux. The zero-dimensional loop model, 
called EBTEL, yields time-dependent solutions for density and temperature in a loop sub- 
jected to additional heating. 

Having no spatial resolution, EBTEL cannot capture the kind of retraction-driven energy 
release we have so far described. Instead it uses a prescribed, uniform volumetric heating 
function, Q{t), to introduce energy to the loop's plasma. We use this to represent retrac- 
tion energy release by specifying a heating profile whose time integral matches SWt/S^. 
Distributed uniformly over a loop 2 Lmm = 120 Mm long (excess accounting for the arcade 
legs) of strength B = Bg = 86 G, the retraction-produced energy amounts to a density of 
117 erg/cm . This is introduced the EBTEL model with a triangular profile over a 50 second 
interval. (Changing the assumed duration has little appreciable effect on the solution, which 
mostly responds to the net energy input). 

Figure [11] shows the evolution of both the DEFT solution (left) and the EBTEL solution 
(right) over their respective simulation times. While the retraction phase, modeled by DEFT, 
finishes in 10 seconds, the evaporation and cooling phase persists for over half an hour. The 
ad hoc heating phase (dashed curves) drives the EBTEL temperature to T ~ 24 MK without 
appreciably affecting the density. Thermal conduction transports this energy rapidly to the 
chromosphere where its is partly radiated (the chromosphere is assumed to be a far more 
efficient radiator than the corona) and the remainder drives evaporative fiow, returning the 
energy through enthalpy fiux. This conductive evaporative phase lasts 600 sec until the 
coronal density reaches rig = 2 x 10^*^ cm~'^ and temperature has fallen to T = 15 MK. 
Thereafter the loop cools through coronal radiation and draining fiow (downward enthalpy 
fiux). 
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Fig. 11. — Two stages in the evolution of the reconnection episode depicted in Fig. |9l The 
left column is the DEFT run covering the retraction through the current sheet. The other 
column is the EBTEL solution for evaporation and radiation inside the stationary flux tube. 
The temperature (bottom), density (middle) and EM (top, actually 6EM/6^) are plotted 
versus time in seconds. The dashed curve in the lower left panel shows the peak temperature, 
while the solid shows the i?M-weighted value. The graph on the lower right shows the DEM 
for the combination of both phases. The grey region shows the DEM of the LTS, from DEFT. 
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The two flare phases can be combined into a single differential emission measure [DEM) 
per unit flux change 

1 r 5EM 



DEM^iT') 



AT 



(5$ 



e[T(t)-T'-iAr]e[r'-iAT-r(t)]rft . (2^ 



where the Heaviside functions, Q{x), restrict the integral to times the temperature falls 
within |AT of the value T'. While it is shown below that this DEM must be interpreted 
with some care, it is plotted in the lower right of Fig. [H] against the same temperature axis 
as the time-profiles. It has a large peak about (T) = 15 MK from the cooling phase and a 
much smaller peak at {T)^^ = 30 MK from the DEFT solutions — the LTS. The cooling loops 
(EBTEL) have a thousand times more EM, per unit flux, than the LTS: 7 x 10^^ cm~^ sec/Mx 
vs. 7 X 10^^ cm~^ sec/Mx. This is due to their much longer life times, and in spite of their 
much lower density (DEFT reaches a density Ue = 8x 10^° cm^^ four times higher than EBTEL). 
Nevertheless, the DEM shows two distinct peaks corresponding to the two phases of post- 
reconnection evolution. These, we propose, are the two distinct temperatures detected in 
flares with SH-LTSs. 

To derive observable properties we must combine the single loop simulation into a com- 
posite. If the assumed discrepant flux A\E' = 3 x 10^^ Mx were reconnected over a typical 
four-minute impulsive phase then the mean flux transfer rate would be $ ^ lO^^Mx/sec. 
Such a value, equivalent to 100 GV , is characteristic of values found from the rate of flare rib- 



bon m otion across magnetograms (jQiu et al.ll2004J : iQiu &: Yurchyshynl 120051 : iLongcope et al. 

2010l ). A single DEFT simulation evolves on far shorter time scales (~ 10 sec, on the left 

column of Fig. [11]), so the SH-LTS would consist of a superposition of individual loop-top 

plugs. If each reconnection patch transferred 5$ = 10^^ Mx, typical of a post-flare loop in 

EUV, the impulsive phase would consist of ~ 300 transfers, occurring at a mean rate of 

one per second. Since each loop-top plug remains visible for ~ 10 sec, roughly 10 would be 

visible at any time, to compose the observed LTS. The emission measure of this composite 

would be 

(^EM\ 



(EM), 



$ 



10*^ cm" 



(29) 



characteristic, albe i t at the low end, of the range observed in SH-LTSs (JLin et al.l Il981 



Nitta fc Yajil 119971 : iLongcope et al.l l2010l : ICaspi fc LinI 120101 ) . It is evident from the top- 
left panel of Fig. [TTl that the EM is very high even at full retraction (i.e. the end of 
the DEFT simultion). It is therefore possible tha t EM persists for some time beyond this 
(iLongcope et al.ll2010l : iGuidoni &: Longcopdl201ll ). If this were to happen the observed EM 
would be larger than our calculation (halted at full retraction) by as much as a factor of two. 

The cooling phase occurs over time scales ~ 1000 sec (right column of Fig. [TTl) . much 
longer than the four- minute flux transfer. In fact, the emission measure of a single loop 
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peaks at 6EM/S^rnx = 4 x 10^^ cm~^/Mx, ten minutes after the energy input. The total 
emission measure from this phase is not therefore found from Eq. fj20|) which apphes only to 
rapid evolution. Instead the net emission measure of the flare as a whole resembles that of 
the individual curves, peaking at 

m&x(EM^^^) = A^ f-?^) - 10^°cm~^ . (30) 

\ / mx 

GOES would register this cooler plasma, with an -EM-wei ghted temperature o f 15 MK, as a 



flux Fi-g ~ 10 ^W/m^ in its low-energy (l-SA) channel ( Thomas et al.l 119851 ): it would be 
an X class flare. 

We have thus concluded that, at least for this model flare, the EM of the SXR compo- 
nent is proportional to A\E', while that of the HXR componen t is proportional to its time- 



deriv ative. This is basically a re-statement of the Neupert effect ( lNeupertlll968l : lDennis fc Zarro 



19931). 



If the evaporative phase persists far longer than the flux transfer, as we have assumed, 
the DEM observed at any instant will not match the time-integrated curve in the lower right 
of Fig. [TT] The SH-LTS will consist of a super-position of loops in all phases of evolution, so 
its DEM will resemble the shaded portion of the plot. The cooling phase, however, evolves 
much more slowly, so its DEM will be more sharply peaked about the temperature presently 
found in all the loops; this would differ from the plotted curve, but still result in a distinctly 
bimodal DEM. The integral under the lower peak would correspond to Eq. fl30|) and will 
be smaller than the integral under the plotted curve, which matches Eq. ( l20l) . The ratio of 
areas will be roughly 1:100, as indicated above. Accounting for the post-retraction phase of 
the LTS could boost the ratio to 1:50 in this particular example. 

The foregoing model X-fiare arose directly from a release of 3 x 10'^^ erg of magnetic 
energy through post-reconnection retraction capable of thermalizing 7 x 10'^^ erg. A LTS of 
(T) = 30 MK would be observable above these cooling loops, although its EM would be 
significantly smaller. The inferred density of the SH-LTS would be ~ 8 x 10^°cm~^, ten 
times greater than the pre-reconnection density. 



7. Discussion 

W e have her e in att empted to show that compression by slow mode shocks, first predicted 



in the iPetscheki ( 119641 ) model of fast reconnection, is capable of producing high-density, 
super-hot loop-top emission like that seen in flares. The density enhancement achieved 
in these shocks is considerably higher than generally assumed on the basis of steady shock 
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models (Rankine-Hugoniot relations). This is due to the coohng of the post-shock plasma by 
thermal conduction during the initial (transient) phase of the post-reconnection retraction. 
This effect has been explored using solutions of one-dimensional shock tube equations and 
two-dimensional thin flux tube equations. 



The simplified shock-tube model led us to an estimate of density enhancement, Eq. ([T| 
in terms of properties of the current sheet and its plasma. Figure [12] uses this relation to show 
the density and temperature enhancements vs. the strength and angle of the fields separated 
by the current sheet. Density enhancements as large as 100 appear plausible for the strongest 
fields [B > 500 G) and largest angles {A6 > 120°). This conforms to observational ev idence 
that SH-LTSs occur only when field strengths are extremely high (jCaspi &: Linll2010l ). 
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Fig. 12. — Estimates of density enhancement over a range of flare parameters. For ambient 
density Ueo = 3 x 10^ cm~^ and temperature Tq = 2 Mk, a range of fields strengths (abscissa) 
and reconnection angles are used evaluate Eq. (IT8l) . (bottom curve). The temperature at 
time of this peak value is plotted in the top panel. The result from simulation in Fig. [7] 
{zr = h/2) is shown with a triangle. The lower triangle is the observed enhancement and 
the upper triangle (on the curve) is the adiabatic value. 



The shock tube model is a natural abstraction of the thin flux tube model for dynamical 
evolution following patchy reconnection. These simplified, time-dependent dynamics turn 



-37- 



out to capture the behavior of even classical steady-state MHD reconnection models, provided 
the reconnecting fields are not perfectly anti-parallel. In all cases shocks, and corresponding 
density enhancements, are direct consequences of energy release. Magnetic energy is released 
by shortening field lines, which naturally compresses the plasma they link. Since the field 
lines shorten at Alfvenic speed the plasma is compressed supersonically (/5 <C 1) leading to 
very strong shocks. 

There is quantitative agreement on shock-produced density enhancement between ana- 
l ytic treatraents which are steady-state (ISowardlll982l : IVrsnak fc Skenderll2005l ) and transient 
(JLin &: Leelll994j : iLongcope et al.ll2009[ ). This agreement occurs because even in steady-state 
models, the response of any single field line to an instantaneous topology change is neces- 
sarily transient. All such models have been based on Rankine-Hugoniot jump conditions 
across the shocks. We have here shown these steady-shock Rankine-Hugoniot relations are 
inapplicable to the transient model since a heat front of asymptotic form is never achieved. 
We fully expect the Rankine-Hugoniot relations will fail in traditional steady-state MHD 
models as well, since those basically capture the same transient dynamics of each field line. 

While analytic treatments might suffer from erroneous assumptions about the shocks, 
we expect time-dependent MHD solutions to reveal the large density enhancements here 
predicted. There have been several two-dimensional simulations of reconnection across a 
current sheet which i ncluded temperature-dependent, fi e ld- aligned thermal conduction nec- 
essary to our model (JYokoyama &: Shibatal 119971 . Il998l : IChen et al.lll999al Jbl). While none 
included a guide field component, the reconnection outfiows consisted of heat fronts ahead 
of isothermal slow shocks — the same structure predicted by the shock tube model. The 
simulat ions used fa i rly larg e initial pressures (/3o ~ 0.1) and thus achieved rather modest 
shocks: I Chen et al.l (Il999bl ) observed a fifteen-fold pressure jump, which would be produced 
by a shock tube infiow at Mi = 2.40. That simulation exhibited a five-fold density enhance- 
ment across the "isothermal shock" , twice the maximum achievable by a switch-off sock in 
steady state but exactly matching shoc k tube pred i ction E g. (IT8|) . Since the anti-parallel 
reconnection {A6 = 180°) simulated by IChen et al.l ( ll999bl ) falls outside the validity of the 
shock-tube model, this agreement must be regarded as fortuitous. 

A more direct treatment of the full problem, including a guide-field, using the full set 
of MHD equations could bypass the limitations of the thin fiux tube model. This model 
requires that the plasma-/3 be small both before and after reconnection. For negligible 
i nitial values (/3n <^ 1 ), the post shock value approaches the limit, (^2 — (16/3) sin^(A6'/4), 
(ILongcope et al.l 120091 ) which exceeds 0.75 for A6 > 150°. Thus reconnection at more acute 
angles, where Mj could be large, cannot be treated using our thin fiux tube model. Indeed, 
there is observational evidence (ICaspi &: LinI |2010| ) that many SH-LTSs achieve a plasma 
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pressure limited by the magnetic pressure {P2 — !)• We therefore hope in the future to use 
full MHD simulations to explore this reconnection regime. 

In both steady or transient scenarios, the total emission measure of the hot reconnection 
outflow will be proportional to the mean reconnection rate $ according to Eq. (I20l) . In 
traditional steady-state models this is equivale nt to an electric field along the X-line (jPetschek 



1964 : 1 VasyliunasI 1 1 9 75l : ISoward &: Priest 



19821 ). directl y related to reconnection microphysics 



This contrasts with patchy reconnection scenarios (JKlimchukI Il996l : iMcKenzie fc Hudson 
I999I ) where the local electric field determines the time required to reconnect each patch, 
but $ is proportional to the rate patches are produced, an unrelated quantity. While it is 
important to resolve this ambiguity if we are to fully understand solar flares, both scenarios 
lead to the same relation between $ and the flare's EM for physical reasons outlined above. 
Moreover, the meai i reconnection rate $ can be directly measured using ob servation of 
flare ribbon motion (JForbes fc Priestlll984j : iFletcher et al.ll200ll : iQiu et al.ll2002l ). The values 
typically obtained are consistent with the EM of typical SH-LTSs. Such measurements have 
been interpreted as X-l ine electric fields by applying the two-dimensional steady state model 
(jPoletto fc Kopplll986l ). Doing so neglects the complexities of flare ribbons and magnetic 



fields , which seem to sugges t a less ord ered, more patchy reconnection process (IFletcher et al. 
2004J : iLongcope et"aDl2007l : [Qiull2009h . 



The constant of proportionality, 5EM/6^, can be derived from global, observable quan- 
tities. It consists of a factor, given in Eq. fl^Tl) . depending on pre-flare conditions and the 
vertical extent of the pre-flare current sheet, multiplied by a dimensionless factor depending 
primarily on the angle A6 between reconnecting field lines which determines the inflow Mach 
number Mi according to Eg. (fTO|). This simple dependence can be illustrated using the 40 
MK SH-LTS studied bv. llongcope et all f l201ol ). which had a constant 6EM/5^ ~ 3 x lO^^: 
it had an emission measure EM ~ 3 x 10^^ cm~^ during reconnection at a mean rate $ ^ 10^^ 
Mx/sec, measured using ribbon motion. Its thirteen-fold temperature increase, from T = 3 
MK to T = 40 MK, requires a shock of Mj ~ 4.7, as would occur for reconnection be- 
tween fields differing by A6 ~ 90° (roughly consistent with the pre-flare magnetic model of 
Longcope et al.ll2010l ). Reading from Fig. [HI this should produce an EM ~ 0.3 times the 
factor in Eq. (12 ip . which must therefore be 10^° secM x"^ cm~^. A current s heet h = 30 Mm 
high with guide field Bg = 200 G (similar to those of iLongcope et al.ll2010l ) would require a 



pre-reconnection density (density inside the current sheet on the un-reconnected field lines) 
rz^n = 2 X 10^°cm~^, t o produce the observed LTS. This is lower than the value quoted by 



Longcope et al.l (120101 ) due primarily to density enhancement above the Rankine-Hugoniot 
value. According to Fig. [Slthe loop-top density in the SH source would be Ue 
five times greater than the maximum permitted by Rankine-Hugoniot relations. 



3xl0"cm-3, 
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Perhaps the greatest puzzle raised by the coUisional model explored here is how it might 
relate to flares (or flare phases) with signiflcant non-thermal populatio ns. Super-hot sources 
ofteri appear after the t e ll-tale power-law HXR spectrum has vanished ( I Alexander fc Met call 
19971 : 1 Jiang et al.l l2006l : iLongcope et al.ll2010[ ). so the collisional and collisionless processes 
may occur separately. It still seems reasonable to assume the same underlying mechanism 
is responsible for energy transfer in both cases. Here we have shown how reconnection can 
thermalize a signiflcant fraction of stored magnetic energy through slow magnetosonic shocks, 
assuming sufficiently high collision rates. Could a related process, at lower densities, produce 
the non-thermal populations observed at other times and in other flares? Except in cases 
of anti-parallel reconnection the slow shock i s essentially a parallel shock, whose collisionless 
mani f estation is st i ll poorly understood (see Parkerlll961t ISagdeevlll966l : IStone fc Tsurutani 
19851 : iQuestI Il988l : iKhabibrakhmanov et al.l Il993l . for overviews of the topic and different 
possible approaches). This seems to be an avenue worth pursuing given that the present 
model outlines a complete energetic chain from current sheet to post-flare loops. 

In any event, the large densities within the SH-LTS are sufficient to thermalize the 
electron population, thereby explaining the absence of a non-thermal population when they 
are present. A troubling discrepancy occurs in trying to explain how bulk kinetic energy, 
produced at the rotational discontinuities and carried almost entirely by ions, can be trans- 
formed to electron thermal energy. Even at the high densities observed the classical rate o f 
ion-electron collision is too low for this transfer to be effective (ILongcope fc Bradshawll2010l ). 
Since all measurements available are of electron temperatures in flares, they must be heated 
somehow, in spite of this theoretical hurdle. This would not be the flrst instance where 
actual collision rates greatly exceeded those of classical Coulomb interactions. Indeed, the 
resolution of this paradox may also provide a clue to the production of non-thermal electrons 
in cases of lower density. 
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